Heavy bottom squark mass in 
the light gluino and light bottom squark scenario 1 

Zumin Luo 2 and Jonathan L. Rosner 3 
Enrico Fermi Institute and Department of Physics 
University of Chicago, 5640 S. Ellis Avenue, Chicago, IL 60637 

Restrictive upper bounds on the heavy bottom squark mass when the 
O ■ gluino and one bottom squark are both light are based on the predicted 

^ ■ reduction of Rb (the fraction of Z hadronic decays to bb pairs) in such 

' a scenario. These bounds are found to be relaxed by the process Z — > 

bbg/bbg, which may partially compensate for the reduction of Rb. The 
relaxation of bounds on the top squark and the scale-dependence of the 
^ ; strong coupling constant are also discussed. 
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As a promising candidate for physics beyond the Standard Model (SM), the min- 
imal supersymmetric standard model (MSSM) has been intensively studied in the 
past years. Recently there has been some interest [1, 2, 3, 4, 5, 7, 8, 9, 10] in a MSSM 
scenario with a light gluino and a light bottom squark. It was first proposed by Berger 
et al. [1] to explain the excess of the cross section for bottom quark production at 
hadron colliders. In this scenario, the gluino has mass m g = 12-16 GeV and the 
| light bottom squark has mass m~ h = 2 - 5.5 GeV. The other bottom squark b , an 

orthogonal mixture of b^ and b r to the light bottom squark, is assumed to be heavy, 
and so are all other supersymmetric (SUSY) particles. We follow the convention in 
Ref. [1] to define 

£ : (~ b ) = ( cos ^ ^o- b \(b R \ m 

[b> ) \ -sin^ cos ^ )\b L ) ■ {1> 

In this scenario the Z decay width Y z is reduced by SUSY-QCD corrections to the Zbb 
vertex [3]. However, some new decay processes contribute positively to Y z and may 
compensate for the reduction. Apart from Z — > gg, which was studied in Ref. [10] 
and found to have a decay width of order 0.1 MeV, there exist two more important 
decays, Z — > bbg/bbg and Z qqgg. The first process is ~ aa s at the tree level and 
has a decay width of 1.9 - 5.9 MeV depending on the sign of sin 26^ and the mass of 
the gluino [7]. The second process is ~ aa 2 s and its decay width is calculated in a 
model-independent way to be 0.75 - 0.21 MeV for mg = 12-16 GeV [6]. Additional 
"sbottom splitting" diagrams [9] raise Y{Z — > qqgg) by about 0.01 MeV. 

By studying the SUSY-QCD effects on R b , Cao et al. [3] showed that the mass 
of the heavy bottom squark should not exceed 195 GeV at the 3a level. Taking 
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into account additional electroweak corrections to the gauge boson propagators and 
including other relevant Z-pole observables and the VF-boson mass, Cho [4] found 
that the mass of b' is further constrained (< 180 GeV) and that one of the top 
squarks (t) should be lighter than 98 GeV at the 5<r level. However, both of these 
analyses ignored the impact of the new Z decays on the electroweak observables. 
Most notably, it was shown in Ref. [6] that Z — > bbg/bbg may affect the measurement 
of Rf,. Indeed, since the gluino is assumed to decay promptly to a bottom quark and 
a bottom squark, the final state in Z — > bbg/bbg will be an energetic bottom jet back- 
to-back with a "fat" jet consisting of a bottom quark and two light bottom squarks. 
The energetic bottom jet can add to Rb if the analysis consists merely of counting 
b candidates. Moreover, if the light bottom squark cannot escape the detector, the 
"fat" jet can also be tagged as a b jet and contribute to Rb [6, 11]. However, since the 
branching ratio of Z — > bbg/bbg is only of order 10~ 3 , current LEP data are not likely 
to be sensitive to it [11] without a dedicated search. Despite this, the Z — > bbg/bbg 
decay plays a nonnegligible role in constraining the heavy bottom squark mass. In 
this Letter we argue that the upper bound on will be relaxed if this process can 
contribute a certain amount to the measured YiZ — > bb). 

The proposed scenario implicitly requires that i?-parity be violated in bottom 
squark decays [4]. Assuming pure lepton-number violation in top squark decays 
(B(t — > rb) = 1), the CDF Collaboration [12] has recently set a 95% confidence 
level lower limit at 122 GeV on the light top squark mass in the framework of 
i?-parity violating MSSM. We will show that the upper bound on can be raised 
above 122 GeV. However, we point out that the CDF lower limit is somewhat irrel- 
evant to the proposed scenario. The top squark decay in the light gluino and light 
bottom squark scenario is not necessarily .R-parity violating or purely lepton-number- 
violating. An .R-parity- violating decay invoked for the bottom squark does not tell us 
anything about how the top squark decays since the two processes arise from different 
terms in the Lagrangian. 

The Zbb vertex in the Standard Model can be written as 

V*(ZW) = -i—L— \r(g b v - g b A i 5 )} , (2) 
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where g\ = — 1/2 + 2/3 sin 2 8w, 9 a = ~ 1/2, and 9w is the weak angle. The SUSY- 
QCD corrections to g\ and g\ have been calculated in the MSSM by a number of 
authors [3, 5, 13, 14, 15]. Using the analytical expressions in Ref. [14], we reproduce 
the numerical results in Ref. [3] . However, we cannot reproduce the numerical results 
in Ref. [5], which also disagree with the other references analytically. 

For numerical evaluation, we take a s (Mz) = 0.120, sin 2 0w = 0.2311, = 4.1 
GeV, m b = 5.5 GeV, m~ g = 16 GeV, Rf M = 0.21569 and r| M (had) = 1.7429 GeV [16]. 
The greatest possible values are chosen for and rrig so that the upper limit on my 
is the least restrictive [3]. The change in R b with respect to its SM value is plotted in 
Fig. 1. For now we do not consider the possible effects of the new Z decay processes. 
The bottom squark mixing angle Q~ h is chosen to satisfy | sin 9- b \ ~ sin 2 0^/3 ~ 0.39 
so that the Zbb coupling is suppressed [1]. We allow | sin 6*^| to vary between 0.30 and 
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Table I: Upper bounds on the heavy bottom squark mass at the 3a level, x is the 
contribution (in units of MeV) of the decays Z — > bbg/bbg to measured Y[Z — > 66) as 
a result of mistagging of the "fat" jets. 

x Upper bound on my (GeV) 
~~ 1 215 
1.5 275 

2 365 
2.5 490 

3 655 



0.45. The partial decay width T(Z -»• 66) is about (1.0, 0.6) MeV for |sin0 s | = 
(0.30, 0.45). Fig. 1 shows that the heavy bottom squark mass my must be less than 
about 200 GeV for R b to be within the 3a bound (SR b > -0.0013) set by the LEP 
experimental value [16]. 

Now suppose that a portion of the Z — > bbg/bbg events cannot be distinguished 
from the Z — > 66 ones and hence contribute x MeV to T(Z — > 66). We shall take 
1.0 < x < 3.0. Note that the decay width for Z — > bbg/bbg is approximately 4.4 MeV 
for = 3 GeV, = 16 GeV and sin^ ~ 0.39 [7]. Since the two bottom quarks in 
the final states will hadronize in the detector, the hadronic decay width T^(had) will 
also be affected. Furthermore, the final-state bottom squark can also contribute to the 
hadronic decay width via i?-parity- violating decays to light quarks [1]. (Possibilities 
of a long-lived 6 and lepton-number-violating decays of 6 to hard leptons are both 
disfavored; see Refs. [1, 17].) Considering also other decays such as Z — > 66, Z — > gg 
and Z — > qqgg, we estimate that the total change in T^(had) is 4.0 MeV+ST(Z — > 66). 
Here 5T(Z — > 66) denotes the change in YiZ — > 66) due to SUSY-QCD effects, as 
discussed above. The upper bound on my can then be relaxed if the additional 
decays contribute positively to R b , i.e., [T SM (Z -> 66) + x]/[r| M (had) + 4.0] > Rf M , 
or equivalently x/4.0 > Rf M . For example, the new upper bound on my is around 
275 GeV at the 3a level for x = 1.5; see Fig. 2, where we take into account the 
additional Z decays and plot 5R b as a function of my for = 5.5 GeV, m~ g = 16 
GeV, sin^ = 0.30 and x — 1.0 - 3.0. The new upper bounds at the 3a level are 
presented in Table I for various values of x. 

SUSY-QCD corrections to g h v and g\ also affect some other related Z-pole elec- 
troweak observables. These are Tz, T(had), Re, R c , and A b , where we use the 

notations in Ref. [16]. Some or all of these observables can also be affected by the 
additional Z decay processes in the proposed scenario. We show in Table II the com- 
parison of the predictions in the SM and in the proposed scenario for my = 180, 275 
GeV, with or without the additional decays taken into account. (Little can be known 
of the effect of those decays on Aj$ and A b .) One can see that raising my from 180 
GeV to 275 GeV does not make the predictions agree much worse with experiments, 
especially after the additional decays are considered. Thus, the constraint on R b plays 
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Figure 1: Change in Rb due to SUSY-QCD effects alone in the light gluino and 
light bottom squark scenario with respect to its Standard Model value as a function 
of the heavy bottom squark mass. Here we take = 5.5 GeV, rrig = 16 GeV. 
The horizontal dotted lines correspond to 2a and 3a bounds set by the LEP R b 
measurements [16]. 
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Figure 2: Change in Rj, in the light gluino and light bottom squark scenario with 
respect to its Standard Model value as a function of the heavy bottom squark mass. 
SUSY-QCD effects and possible impact of the decay process Z — > bbg/bbg are com- 
bined, x is the contribution (in units of MeV) of the process to measured Y(Z — > bb) 
as a result of mistagging of the "fat" jets. From bottom to top, x takes the values 
1.0, 1.5, 2.0, 2.5 and 3.0. Here we use m~ b = 5.5 GeV, m~ g = 16 GeV and sin^ = 0.30. 
The horizontal dotted lines correspond to 2a and 3a bounds set by the LEP Rt> 
measurements [16]. 
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Table II: Deviations of the predictions for some Z-po\e observables from experimental 
values [16]. MSSM 180j 275 (1) [(2)] denotes the prediction in the light gluino and light 
bottom squark scenario for my = 180, 275 GeV before [after] taking into account the 
additional Z decays. Here we use m~ h = 5.5 GeV, m g = 16 GeV, sin^ = 0.30. 



Observable SM MSSM 180 (1) MSSM 180 (2) MSSM 275 (1) MSSM 275 ( 2) 





0.61(7 


-0.52a 


1.87(7 


-1.00a 


1.39a 


r(had) 


-0.75a 


-2.05(7 


-0.05(7 


-2.60a 


-0.60a 


Re 


-1.20(7 


-1.82(7 


-0.87(7 


-2.08a 


-1.12a 


R c 


-0.19(7 


-0.11(7 


-0.23(7 


-0.07a 


-0.20a 


4(0,6) 

^FB 


3.18(7 


3.39(7 




3.43a 




A b 


0.69(7 


0.85(7 




0.88a 





the most important role in setting an upper bound on my, which we have already 

shown can be relaxed to about 275 GeV if the process Z — > bbg/bbg contributes 1.5 
MeV (about one third of its own decay width) to the measured T(Z — > 66). With 
my ~ 275 GeV, we argue below that the upper bound on mi can also be relaxed. 

The constraint on m ( - originates primarily from SU(2)l gauge symmetry and the 
constraint on oblique corrections, i.e., electroweak corrections to gauge boson prop- 
agators. These corrections are parametrized by ASz, ATz and AM\y in Ref. [15]. 
It has been shown [4, 15] that the dominant contributions to these parameters come 
from the bottom and top squark corrections (denoted by AT) to the T-parameter 
[18]. Since T is defined in terms of vacuum polarization amplitudes of the SU(2) L 
gauge bosons, only the left-handed components of the squarks can contribute to AT. 
For given bottom squark masses rag, my and left-right (L-R) mixing angle 9~ b) large 
L-R mixing for the top squarks is favored so that the left-handed component of the 
lighter top squark mass eigenstate is relatively suppressed and therefore a small AT 
can be obtained [4]. The L-R mixing for the top squarks is parametrized by A* ff in 
Ref. [4], where it is constrained to be no less than 300 GeV for my = 180 GeV in 
the light gluino and light bottom squark scenario. This leads to a very restrictive 
upper bound (98 GeV) on raf For A l cS = 300 GeV and a consistent set of masses 
and mixing angles (rag = 5.5 GeV, my = 180 GeV, mi = 98 GeV, rap = 340 GeV, 
9- b = arcsin0.30 ~ 17.5°, 0~ t = 49.2°), we find AT ~ 0.081. Larger values for AT are 
obtained if mi > 98 GeV. Without going through a similar \ 2 fit procedure to that 
in Ref. [4], we regard 0.08 as a threshold value for AT below which oblique correc- 
tions are assumed to be acceptably small. For the same A l eE but with my raised to 
275 GeV, we find that smaller values for AT can be obtained for mi as large as 210 
GeV; see Table III. This implies that the upper bound on m t ~ may be considerably 
relaxed if the process Z — > bbg/bbg can contribute a certain amount to the measured 
F(Z -> 66). 

One can even tolerate much larger values of AT if one is prepared to consider 
heavier Higgs boson masses and small positive changes in S in precision electroweak 
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Table III: Examples for m- t > 98 GeV and AT < 0.08. We take m~ b = 5.5 GeV, 
my = 275 GeV, sin% = 0.30. 



m~ t (GeV) 


m~ t , (GeV) 


h 


AT 


100 


357 


31.9° 


0.061 


165 


368 


38.1° 


0.070 


210 


387 


46.4° 


0.079 



fits [19, 20]. In one example [20], by letting the Higgs boson mass rise to about 1 TeV, 
one can accommodate AT as large as 0.5 at the price of AS = 0.3 while retaining an 
acceptable fit to the data. 

Finally we comment briefly on the scale-dependence ("running") of the strong 
coupling constant a s in the light gluino and light bottom squark scenario. In the 
context of this scenario, we showed in a previous work [8] that a s (Mz) fell in the range 
(0.130 - 0.135) ± (>)0.003 if extrapolated from low mass scales (such as m&) and in the 
range (0.123 - 0.131) ± 0.005 if extracted directly from the Z-pole observable T^(had). 
These two ranges overlap with one another and no clear-cut decision could be made 
in favor of either the Standard Model or the proposed SUSY scenario. However, we 
were only able to evaluate Y[Z — > gg) very roughly and did not take into account 
the decay process Z — > bbg/bbg at that time. With both of these processes now well 
understood and considered, the possible range of the extracted a s (M z ) from T z (had) 
turns out to be different. For numerical purposes, we again use a s (Mz) = 0.120, 
m - b = 5.5 GeV, m- g = 16 GeV, sinflg = 0.30 and take 120 GeV < my < 655 GeV. 
Here the lower bound on my is based on a recent analysis by Berger et al. [21]. The 
SUSY-QCD correction to T(Z -> bb) is then 5T(Z -> bb) = -(1.3 - 5.6) MeV. We 
neglect the relatively small changes in 5T(Z — > bb) due to different choices of a s (Mz)- 
Assuming the contribution of the new decay processes to T^(had) is 4 MeV, the net 
change in the predicted T z (had) is then (—1.6 — 1-2.7) MeV. This change can be 
accounted for by tuning a s (Mz) to lie in the range (0.118 - 0.126) ± 0.005. It is 
lower than the extrapolated range from low mass scales, but not at a statistically 
significant level. The heavy sbottom mass should be better constrained to reduce the 
indeterminacy. 

To summarize, we have shown that it is possible to circumvent restrictive upper 
bounds that have been placed on the heavy bottom squark mass when the gluino and 
one bottom squark are both light. The reduction of Rb predicted in such treatments, 
which would conflict with data, can be compensated by a contribution to the b quark 
production cross section from the process Z — > bbg/bbg. In such a case one expects a 
fast b quark jet to be accompanied by a "fat jet" associated with the hadronization 
of the bg or bg system. The relaxation of the bounds on the heavy bottom squark 
mass is accompanied by a corresponding relaxation of bounds on the top squark mass. 
The situation of the scale-dependence of the strong coupling constant a s becomes less 
favorable to the light gluino and light bottom squark scenario after the decay process 
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Z — > bbg/bbg is taken into account. However, improved understanding of the effects 
of the new decay processes in the scenario on Rf, and r^(had) is needed for a clear-cut 
conclusion. 
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